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^ . Quantum Mechanics places limits on achievable transverse beam spot sizes 

of particle accelerators. We estimate this limit for a linear collider to be 
Ax > hcf /E5q where / is the final focal length, E the beam energy, and 
5q the intrinsic transverse Gaussian width of the electron wave-function. <5 is 



determined in the phase space damping rings, and we find 8 ~ yhc/eB where 
B is the typical wiggler magnetic field strength in this system. For the NLC 
5q ~ 25nm, and Ax ~ (9(0.06) nm, about two orders of magnitude smaller 
than the design goal. We can recover a crude estimate of the classical result 
when we include radiative relaxation effects. We also consider a synchrotron 



and obtain Ax > yhcf /E ~ O(1.0) nm, We discuss formulation of quantum 
beam optics relevant to these issues. 



1 Introduction 



Particle accelerators are designed and built, based essentially upon the classical theory of 
point charges interacting with electromagnetism. Nevertheless, particles are described by 
wave-functions, and diffractive limits must exist as to how well they can be localized in 
a given optical apparatus. The first quantum mechanical effects to arise in a potentially 
limiting way might be expected to be diffractive in nature. In this paper we take a first 
look at the problem of estimating the quantum diffractive limits of accelerators. We 
begin with an important system, an NLC-class machine. We are inspired to consider 
this because the desired goals for the NLC beam spot size are ambitious. To achieve the 
desired luminosity requires a ~ 5 nm beam spot in one transverse dimension, (the vertical 
or y direction in the NLC reports 0). We will find that this criterion is about two orders 
of magnitude above the quantum limit. Indeed, we will describe how to estimate the 
classical design result for the beam spot size itself from quantum mechanics, obtaining 
rough agreement with the NLC specifications. 

We obtain a conceptually simple result. The diffractive limit on the beam spot size in 
the x± direction is given by the Rayleigh formula for a (massless) wave of energy E which 
has passed through an effective "aperature" 5o and focused over a focal length ~ /. That 
is: 

Az>M (1.1) 
E6q 

where / is the final focal length, E the beam energy (the result varies somewhat in a 
compound lens system, see Section 3). We emphasize that the "aperature" 5$ is not a 
mechanical aperature, e.g., it is not the beam pipe size. 5$ is actually the initial state 
Gaussian width of the transverse quantum wave-function as it enters the linac upstream 
from the damping rings, where the wave-function has been prepared (we assume an "ideal 
linac," in which there is negligible further synchrotron radiation downstream; this is not 
necessarily a good approximation, and corrections to the effective S are expected). The 
initial state can be considered to be an ensemble of particles, each in simple harmonic 
oscillator (SHO) transverse wave-functions, where the Gaussian envelope (groundstate) 
width is determined by the damping ring wiggler system. This is given by: 




(1.2) 

where B is the typical magnetic field in the damping system, of order 1 Tesla. Taking 
/ = 2m, E = 250 GeV, B = 1 Tesla, yields Sq ~ 25.7 nm, and thus, Ax > 0.062 nm as 
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a diffracive limit. Hence, the NLC would appear to be safely above the quantum limit 
by about two orders of magnitude. We remark, however, that this is the extremal lower 
limit which saturates the Heisenberg uncertainty relationship, and holds in our idealized 
limit. 

In general the individual particle initial state is an excited SHO transverse wave- 
function of average principle quantum number n. This increases the expected diffractive 
spot size to Ax ~ \ffif / 'E5q. In fact, since n oc 1/H it is easily seen that this result is 
independent of h, and therefore should be equivalent to a classical derivation of the beam 
spot size, n can be crudely estimated from radiation relaxation following the original 
arguments of Sands |§ , and others. This yields a result of Ax ~ 2 nm, roughly consistent 
with the NLC design report calculations for the vertical beam spot [|I|]. 

The subject of quantum beam dynamics for particle accelerators is fairly novel 0. 
Gaussian optics is a preferred formalism for tackling the problem considered here. Presently 
we construct a transverse Gaussian wave-packet, with a longitudinal plane wave struc- 
ture, and propagate it through an optical system. Gaussians extremalize the Heisenberg 
uncertainty relation, and they are also the groundstate solutions in continuous linear fo- 
cusing channels, and e.g., magnetic lenses, wigglers, to a reasonable approximation, etc., 
and can be approached by synchrotron radiation relaxation ||, 0], ||, ||. Remarkably, 
Gaussian transverse wave-functions, which solve the quantum Schroedinger equation for 
propagation through the optical system, (neglecting synchrotron radiation), are controlled 
entirely by the classical lens matrices of the system. While Gaussian optics is a standard 
formalism in treating electron microscopy 0, ||, ||, to our knowledge, the behavior of 
a quantum Gaussian beam in a synchrotron has not been previously formulated, and we 
will indicate the self-replicating solution to a synchrotron by an application of lens matrix 
methods. 

First consider the problem of a relativistic electron wave-function passing though a 



lens. Spin is an inessential complication [1C], [ II], so we can use the Klein-Gordon (KG) 



equation. Assume for simplicity that there is only one spatial transverse dimension, 
x_i_, and let z be the longitudinal spatial dimension. In the KG equation we include a 
transverse simple harmonic oscillator (SHO) potential term which is dependent upon z 
(For the analysis, we set h = c = 1): 

d 2 <p + m 2 + K{z)x\<f> = (1.3) 

Then, with = exp[— i(Et — p z z)]4> and E 2 = p 2 + m 2 , the KG equation becomes the 
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transverse Schoedinger equation: 

*|# + ^W-^# = (1.4) 

where: 

K = K/E. (1.5) 

This is a standard construction in optics ||, and (f>(z) has the conventional interpretation 
with z replacing time. The parameter K(z) is ^-dependent, corresponding to the finite 
longitudinal structure of the lens system. For a single thin lens we take K(0) = for 
z < —5z and z > 0, and K = Kq for — Sz < z < 0. 

Let us now postulate a Gaussian form for the wave-function centered at the transverse 
position x±, carrying a transverse momentum p±: 

4> = exp (—-A(z)(x — x±) 2 + ip±x + (1.6) 

In this expression, A(z) is the complex Gaussian kernel, x± and p± are real, and C simply 
parameterizes the overall normalization. Hence, the Gaussian wave-function has four real 
parameters. After substituting this anzatz, the Schroedinger equation, eq.fll.4p, yields the 
following equations of motion for the width: 

= (1-7) 
and for x± and p± we obtain the classical Hamilton equations: 

dp± dx ± p ± 

-dz~ = ~ Kxr > ~~dz~ = ~E ' (L8) 
Note that the centroid x± and centroid momentum p± motions are decoupled from that 
of the Gaussian kernel A(z) and vice versa, (anharmonic effects would generally couple 
these quantities). Moreover, the boundary conditions on A(z) and of the centroid x± 
and centroid momentum p± are independent . Note that the last equation is just the u z- 
velocity" expressed in terms of the momentum for a particle of "mass" E. Remarkably, 
eq. (|1.7|) can be seen to be equivalent to the classical Hamilton equations f| by identify- 
ing A(z) = iP(z)/X(z) where P(z) and X{z) are generalized (complex) momenta and 
positions which satisfy the eqs.(|1.8|). 

Now, we impose an initial condition at z = —L that the particle has been prepared 
into a transverse Gaussian wave-packet, specified to have a pure real width 5 given by: 

A = Re[A(—L)} = l/(5 ) 2 ; Im[A(-L)\ = (1.9) 
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We assume that the centroid of the initial wave-packet is moving parallel to the z-axis, 
thus p± = 0, and x± = Xq is initially arbitrary at z = — L. 

The wave-packet enters the lens at z = —5z and exits at z — 0. Upon entry of the 
lens A(— Sz) is given by the free drift solution of eq.( |1.7| ) from z = —L to the lens, over 
the drift distance L: 

where we assume a thin lens, Sz/L << 1. 

In the thin lens, to a good approximation for small 5z, we have from eq. ( |1.7| ): 

A(p)=A(-6z)+i8zK (1.11) 

Here we neglected the term i5zA 2 /E in the differential equation for A(z) which only gives 
negligible free particle spreading in the lens. Moreover, the classical centroid motion of 
the wave-packet is found from eq.( |1.8|) : 

Rl(0) = -K x 5z; zl(0) = x . (1.12) 

Upon exiting the lens the particle propagates again in free space a distance I with K = 0. 
Hence we find: 

and: 

P±(i) = -K x 5z; x ± (£) = x -p ± i/E. (1.14) 

Note that the classical trajectory of the off-axis particle is deflected back toward the lens 
axis, x± = 0. 

The focal length, / is defined such that x±(f) = 0, hence: 

f = J^T (1-15) 

KqOZ 



The kernel of the wave-packet can now be obtained by solving eqs. ( |Q , 1.1 U| , |l . 1 1] , |1 . 1 3| ) 
recursively to obtain: 



A(£) 



l-i[L/{E{5 Q f)-{5 fE/f] 
(5 ) 2 (1 - t/f) + i{£/E + L/E - Li I Ef) 



;i.i6) 



Note that the Gaussian kernel has an imaginary part which changes from positive 
(focusing) to negative (defocusing) upon passage of the geometrical focal length, L > f. 
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Figure 1: The numerical evolution of a Gaussian wave-function with L = 3.0 evolving 
through a lens; E = 1.0, K = 10.0, 5z = 0.1 and A(-3.0) = 3.0, hence 5 = 0.5 77, and 
the geometrical focal length, / = 1. (a) dashed (black) line is x±(z) crossing the axis at 



the geometrical focal point; (b) dotted (red) Gaussian width l/jRe(A(z)) which focuses 
at f q ; Also shown are Re(A(z)) (green solid) and Im(A(z)) (blue dot-dashed); note the 
Im(A(z)) receiving a positive kick upon passing through the lens. 

Thus, the transverse probability distribution becomes: 



In Figure 1 we give a numerical integration of the Schroedinger equation in the pre- 
ceding discussion, which confirms the validity of our solution. Note that for finite L the 
Gaussian width is focused to a minimum at z — f + f 2 /L + .... In the limit L —>■ oo the 
transverse size of the wave-packet reaches a minimum at the focal point £ = f, where the 
new effective transverse size is: 



This is the usual Rayleigh diffractive minimum, fX/a if we regard a ~ 5 as an "effective 
aperture size" through which the beam has passed, and E = p z c = hc/X, the usual 
quantum wavelength of the particle. 

What if the initial prepared wave-function is not the groundstate of a SHO (pure 
Gaussian), but is rather an excited eigenstate of principle quantum number nl Hence, at 
z = —L, neglecting x± and p±, we assume: 





(1.17) 



5\£) = (5 ) 2 (1 - £/f) 2 + (£/E6 ) 2 [l - L{£ - f)/f£f 



(1.18) 




(1.20) 
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where H n (£) is the nth Hermite polynomial. 

First, we note that the Gaussian solution eq.( |1.6| ) contains the generating function for 
Hermite polynomials [ |T2|j : 

= exp --A(z)(x) 2 + tp ± (z)x + C PT " j W(A(z)/2)x±(z)Y 

V A 7 \n=0 n - 

(1.21) 

For a freely drifting particle, if we choose p±(— L) = 0, and initial x±(— L) = xo, in 0, 
then we see that our solution ip(z) is determined for any z: 

^W = ^)k=o (1-22) 

After passing through an arbitrary lens system, the solution for ip(z) becomes messy, and 
in general x±(z) and p±{z) are arbitrary, and we cannot so easily differentiate with respect 
to x to pull out our solution. However, both x±(z) and p±(z) are proportional to x by 
the linearity of the Hamilton equations. At a focal point we have x±(f) = (for any xq, 
owing to linearity), and p±{f) oc xq. Hence, the solution at a focal point simplifies: 

nf) = a^w)l*°=° 



d n ( 1 \ 

exp i--A(z)(x) 2 + ip±(z)x + Cj | a ,,^ n 



dx, 



o 



oc x n exp (-^A{z){x) 2 ^J (1.23) 

and thus, the arbitrary solution is focused to a Gaussian times a power of x. This gives 
a focal spot size: 

Ax = # (1-24) 
Now this result may seem counterintuitive; we are starting with a broader initial distri- 
bution by the factor s/n, and we might guess that this would produce a smaller focal 
point by an amout \j\fn. The wave-function, however, is not smooth in x, i.e., the 
Hermite polynomial yields a distribution of transverse momentum, and the initial state 
has "ears", each of typical Gaussian width 5q, but displaced off the optical axis by yfn . 
These produce the ^Jn enhancement of the focal spot. Yet another way to see this is to 
note that one can make a classical off-axis centroid motion of the groundstate Gaussian 
by superimposing large n states, and the Gaussian width will yield the minimal f/SoE 
result. Of course, the quantum state of interest to us will typically have a large value of 
n determined by radiative relaxation. We consider this in the next section. 
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The actual linear acceleration phase is inconsequential to this result. The above discus- 
sion assumed a uniform drift in the longitudinal ^-direction, i.e., constant energy E. If the 
particle is accelerating linearly, then E becomes z-dependent, E(z) = (Ef — E )z/L + E . 
It is easily seen that the only effect on our solution is to replace L by L hi((Ef — E )/E ), 
where Eq is the initial energy, Ef the final energy, and E in the above expressions is 
everywhere replaced by Ef. For the first NLC, we have Ef ~ 250 GeV >> Eq ~ 2 GeV. 
The linear acceleration phase is thus equivalent to free drift through an effective distance 
of Lln(Ef/Eo) ~ 45 km where L — 10 km. The amount by which a wave-packet of 
initial size of 25 nm spreads throughout the NLC acceleration phase is about a factor of 
6. However, this spreading is irrelevant to computing the final diffractive limit as seen in 
eq. (|1.19|) where the free drift length L has completely cancelled from the expression at the 
classical focal point, and only the initial quantity S (together with the local quantities / 
and E) controls the diffractive limit .[] 

Thus, the ultimate diffractive limit is controlled by the initial boundary conditions on 
the wave-function size, i.e. by 5q, and not by the intervening unitary lens system. What 
in general determines 5q7 For the NLC the initial wave-function, as well as the initial 
classical distribution, is prepared in the "damping rings." Damping rings are essentially 
a system of magnets arranged as wigglers which induce synchrotron radiation and cool 
the classical beam bunches of electrons. They are designed to produce roughly a four 
order of magnitude reduction in one of the transverse dimension phase space volumes, 
i.e., ~ AxAp x (the transverse emittance). As the system cools classically, it is also 
relaxing quantum mechanically. This occurs because the particles in the wiggler chain 
experience a transverse SHO potential, and synchrotron radiation pushes highly excited 
wave-functions toward the Gaussian groundstate in this potential ||. However, there are 
also re-excitation transitions which eventually come into equilibrium, and a typical average 
SHO principle quantum number is established. While this is certainly an oversimplified 
view of the actual system, we will use it as a starting point to estimate Sq. 



1 We remark that the proper way to view the quantum spreading in the transverse phase-space is to 
use Wigner functions, which depend upon both x and a quantum momentum p. The Wigner function 
isocontours deform in a manner that is conformal to the classical emittance envelope, so while the wave- 
functions spread in x the Wigner envelopes actually shear in x and p and remain contained in the 
transverse phase-space. jl5| 
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2 Magnetic Focusing and Damping 



We now sumarize the details of the motion of a transverse wave-packet in a magnetic field. 
This is discussed in detail in the classic work of Sokolov and Ternov [|J. We will use the 
more transparent WKB approximation, expanding about the classical radius of motion. 
Hence, one should use caution in comparing solutions, e.g., principle quantum numbers 
refer to different things. For example, large n in the usual framework Q] corresponds to 
small n, but large classical radius j z presently. 

Consider a particle moving in a planar orbit in a uniform magnetic field, aligned in the 
z direction, B = zBq in a cylindrical coordinate system (r, <fr, z). The vector potential can 
be chosen as = tBq/2 with A r = and A z = 0. We examine the transverse motion 
of a relativistic electron in the plane z = 0. For an anzatz of the form e~ lEt / h i/j(r, (p), the 
KG equation becomes: 



-E 2 + m 2 - 



d 2 



1 d 



+ 



i d 
r d(b 



eA, 



ip = 



(2.25) 



where a possible momentum component p z has been set to zero. Consider a state of with 
a large "pseudo-angular momentum," £, and scale out a factor of 1/ y/r: 



^ = —e ii <i>x(r,t). 



(2.26) 



{£ is not the physical angular momentum because it is gauge dependent due to presence 
of the vector potential; the physical angular momentum in the present case is 2£, as we 
will see below). Hence: 



E? + m?-J*!-+(-i- -erB V + — 
dr 2 \ r 2 / 4r 2 



(2.27) 



This now has the apparent form of a one-dimensional Schroedinger equation with an 
effective potential: 



V(r) 

The potential has a minimum at: 



(2.28) 



{ eB 



1/2 



/ 21 
7b~ 



(2.29) 
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where the latter expression corresponds to I >> 1. Consider henceforth £ >> 1. We 
consider small radial fluctuations around the large orbital radius _R as r = _R + x and 
expand: 

r ft 2 i _ i 

X = (2.30) 



-E 2 + m 2 -^ + V(Ro) + \x 2 V"{R,) 



and: 

V(Rq) = 2eB i V"{Rq) = 2e 2 B 2 (2.31) 

Thus the high orbital angular momentum Landau levels are approximate eigenstates of 
the SHO potential defined by V(x) = \x 2 V"{Rq), or K = e 2 B 2 /2E. The states are 
labelled by (n, £) where n is a principle SHO quantum number; the energy eigenvalues of 
these levels are given by: 

E 2 = m 2 + eB (2£ + n + i) (2.32) 

In the presence of the gauge interaction the physical angular momentum is L z = —id/d<j)+ 
eRoA^Ro) . Hence, the physical angular momentum is: 

j z =£+ -eBoR 2 = 2£ (2.33) 

where we use the explicit solution for Rq from eq.( [2.29 ) in the latter expression. Therefore, 
to make the classical correspondence, we identify the angular momentum with that of an 
entering beam particle of momentum p^, to obtain Rop^ = j = 21. This yields consistency 
with the familiar expression for the classical orbital radius and the total energy: 

Ro = ^r E 2 = m 2 + eB ( 3z + n + h (2.34) 

Transitions that increase n, but decrease j z are allowed; hence synchrotron radiation can 
be excitatory as well as relaxational. The fact that the energy is degenerate, depending 
upon the combination j z + n is a consequence of the symmetry in the choice of the 
classical orbital center. (Note that the solution formed with the anzatz e~ a ^ for large I 
is actually a solution of vanishing physical momentum; it is a zero-mode associated with 
the translational invariance of the center of the particle's orbit). 

The groundstate in the transverse dimension is a Gaussian with A = \eBo\, given by: 

' (eB /hc)- 1/2 . (2.35) 



Re(A) 



For a typical field strenth of 1 Tesla we obtain 5q ~ 25 nm. The "spring constant" is 
0(e 2 ), hence we say that the dipole magnet is weakly focusing (for quadrupoles V"(R ) ~ 
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eB E/a, where a defines a gradient, hence "strong focusing"). This description applies 
to wigglers, even though the dipole magnet field is alternating in z, if the magnitude of 
the B field is roughly constant. 

It has been known for a long time that an equilibrium between deexcitatory and exci- 
tatory transitions for a particle in a damping system (or synchrotron) will be established 
Hj, and there will be an equilibrium value of n. This value is roughly estimated as follows. 
The typical energy of synchrotron radiated photons is H , j|] : 

*- ~ 5 U (2 ' 36) 

A unit step in a quantum number n or j z produces only a small energy change, ~ eB /E ~ 
1/Rq. The dipole approximation selection rules imply large allowed changes in j z , but 
only unit steps in n: 

E { E \ 3 / E \ 3 
Aj z ~ -— — f — ) ~( — ) An~±l (2.37) 



RqcBq \m e J \m 

Bear in mind that we are treating n as the principle quantum number in the WKB focusing 
channel defined by expanding about Rq, and dipole transitions involving the operator A- V 
will change n by a unit (these can be excitatory). Then Aj z /R is essentially the change 
in longitudinal electron momentum, imparted to the photon. In transitions, though Rq 
changes, there is no sudden translation in the transverse position of the electron wave- 
function, only a transition in motion, i.e., the virtual center of the orbit changes @, ||. 

Over a radiative energy loss time interval the number of emitted photons is: 

(2 - 38) 



The principle quantum number n undergoes a random walk by roughly ^/n^, hence the 
equilibrium n is of order ~ y/n^. Using E ~ 2 GeV, and B Q ~ 1 Tesla, whence Ro ~ 6.6 
m, we find n 7 = 1.12 x 10 6 and n ~ 1.06 x 10 3 . 

Hence, our diffractive limit is now increased by \ffi ~ 0.33 x 10 2 , and we thus have 
a beam spot size \ffi x 0.06 ~ 2.0 nm. Why is this result so close to the design goals of 
the NLC that are obtained by classical physics? Indeed, we believe that this result is a 
quantum derivation of the classical result! The quantum number n scales as 1/h, while 
our diffractive limit scales as Ax oc y/K, hence the product ^/nAx is independent of ft. 
This, moreover, assures us that the ultimate quantum limit is of order 1/y/E smaller than 
the minimal classical analysis. (We note that the Oide effect || may be understood as a 
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blowing up of n in intense final focus magnets, where large transverse energy photons are 
radiated). 

A more detailed discussion of synchrotron radiation relaxation is beyond the scope of 
the present paper. Excellent treatments can be found in ||, |J, and the pioneering work 
of §. 

3 Quantum Particle in a Synchrotron 

The solution to the Schroedinger equation for passage of a free particle through a lens, 
eq. flnp , can be completely described by the simple classical optical matrix methods. If 
one passes a classical ray moving in the z direction through a lens system, the outgoing 
state of the transverse x canonical variables may be written as pi: 



The unimodular matrix M. for a compound sequence of lens elements is the corresponding 
sequential product of the individual matrices of the elements. For example, a sequence 
of free propagation (distance L), followed by defocusing lens (focal length — /), followed 
by a space a, followed by a focusing lens (focal length /), followed by free propagation 
(distance t) yields the result: 

aL a£ aLl 



( x{£) \ = ( 1 + J- ! F (L + a + £) + f-f-!f\/ x (-6z) 
{ P{i)/E ) \ -f l -J-f" )[ P(Sz)/E 



(3.40) 



The zero of the (11) matrix element in £ = F = f + f 2 / a implies the system is net focusing 
with composite focal length F (e.g., see ref. JL3|). 

The effect of this particular lens system in quantum mechanics, e.g., on the Gaussian 
kernel A as defined in eq. Ql.6p , can be easily derived from the Schroedinger equation: 

A(f) = A {l-a/f-aL/p)+iEa/p 

U 1 + a/f -a£/f 2 + iA /E[L + a + £ + aL/f -ai/f -(a£L/f 2 )} { ' ' 

The focal length, F, is where Mu — 0, and, at the focal length we obtain the width: 



m = where: R < A (°v = ^ ( 3 - 42 ) 



ii 



This result is the minimal diffract ive quantum limit for the composite lens system, and it 
is again determined by the initial width of the quantum state.0 

Of course, beyond do, there is actually no new information in the above formula for A 
than is already present in the lens matrix for the classical ray optics. If the lens matrix is 
Aitj, then we see by comparison with eq. fl3.41p the general result for the Gaussian kernel 

(0, §): 

Theorem I: A out = — — — — - — — (3.43) 

M11 + iMi 2 A in E 



That is, if we write the out-amplitude as: 

= -iE ( 

V 

then: 

Ma Mu 




A out = -iE (^-) (3.44) 



(3.45) 




^ M 2 i M22 

From this result we can easily derive the quantum limit on beam size at the classical 
focal length. The classical focal length occurs where Mn — 0. At this point we have 
M12 = — M 2 \- Hence we readily find: 

E 2 Re A 1 

Theorem II: Re A* = , j0 r . - - . — — — — = . - . (3.46) 

; Mj 2 [(ReA ) 2 + (ImA ) 2 ] (Ax) 2 v ; 



If A = Ti/5q is pure real, then: 



Mw F 



noting that M\2 = F ~ / is a length scale comparable to the focal length at the clas- 
sical focal length, e.g., F = f 2 /a in our previous compound lens example (see previous 
footnote). 

Now, if the magnet system is periodic, as in a synchrotron, we expect quantum states 
that are approximately periodic solutions in the matrix. Periodic solutions must be eigen- 
states of the matrix M. Consider first the motion within a very thick lens, i.e., a contin- 
uous transverse SHO potential. For an infinitesimal displacement in the ^-direction, the 
lens matrix is: 

Msho=( L 6 ?) (3.48) 



2 Here we might imagine taking a — > 00 holding / fixed to cause S(F) — ► 0; however, for a >> f the 
longitudinal dimension (Az) of the focal point becomes small as ~ f/a owing to the aL/f and alLj f 2 
terms in M12, and the finite longitudinal distribution of the beams becomes problematic; we have not 
looked in detail at optimization of this. 
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which is unimodular to 0((Sz) 2 ). The eigenvalues of Mh are A± = 1 ± i8zJ Kj E. The 
stable, quantum solutions in the lens are therefore the eigenvectors: 



A 



± 



-iE 
A 



Mu M12 
M 2 i M22 



-iE 
A 



hence we find: 



Ac 



±VKE 



(3.49) 



(3.50) 



Hence, the stable solution in the linear focusing channel is, indeed, the Gaussian ground- 
state solution in the SHO potential. We have simply recovered the usual Gaussian ground- 
state in this limit. 

Consider now a "synchrotron," i.e., periodic magnet lens lattice based upon the above 
lens configuration. We assume an infinite series of alternating dipoles with spacing a and 
focal lengths ±/. Replacing L = a — i in the matrix elements Mij of eq. (|3.40|) gives the 
lens matrix for the synchrotron: 



M 



1+2. 



a( 

7 



2a + ^f-2f 



a 2 £ 1 a£ 2 

~F ^ J 1 



(3.51) 



T x / 

The condition that we have a periodic solution is: 

M 2 2A e - iM 2 iE 



7^7 



At 



Mn + iMi 2 A e /E 



(3.52) 



Using det M = 1 we find: 



At 



-iE 



2M 



12 



A^22 - Mu ± ( (M22 + Mn) 2 - 4 



i/2- 



(3.53) 
therefore 

(3.54) 



Stable quantum solutions (solutions that are normalizeable Gaussians for all 
require: 

- 1 < ^Tr(M) < 1 

This is, of course, the familiar stability condition for the classical motion. Note that 
Tr(M) = 2 — a 2 /f 2 , which is I independent, thus when the condition is met for particular 
choices of / and a it holds everywhere. The stability condition is the usual one, / > a/2. 

The solution for A(£) is: 



A(£) 



E 



[(2/ + a)(l -<//)+*»//] 



a" 

W 2 . 



1/2 



1 \ 1 + Tf-1 



(3.55) 
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where f > a/2 and < i < a. 

Consider the special case of a system in which f = a. We see that the minimum 
Gaussian width occurs at £ = a, given by: 

min{l/ReA{£)) = = {Ax) 2 (3.56) 

This implies that the minimum achievable beam spot size in a synchrotron is ~ \J f / E. 
If a focusing magnet of focal length /' is inserted into the synchrotron magnet lattice 
and, then we obtain the minimum spot size ~ f'/EAx ~ f'/y/Ef. There is no initial 
parameter Sq since we have assumed that synchrotron radiation relaxes the quantum state 
into the stable, periodic solution. Here we see a potential advantage of a linear collider 



over a synchrotron, in that the linear collider has a much larger Sq ~ y f / E prepared in 
the low energy damping ring, which makes the quantum diffractive limit smaller, ~ f / ESq, 
while in the synchrotron So ~ \J f / E, where E is the larger beam energy, thus giving a 
larger diffractive limit ~ y/f/E. 

4 Conclusion 

In conclusion, we estimate the minimal quantum beam spot size achievable in a linear 
collider to be given by: 

Ax > M (4.57) 

ESq 

where / is the final focal length, E the beam energy, and So is the initial transverse size 
of the wave-functions prior to acceleration. This may be viewed as a direct transcription 
of the Heisenberg uncertainty principle. Sq is prepared in the synchronous damping rings, 



typically wigglers, and 5o ~ 1/ yeB/hc where B is the magnetic field strength. We have 
for an NLC-class machine, B ~ 1 Tesla, / ~ 2 m, E ~ 250 GeV hence Sq ~ 25 nm, and 
Ax ~ 0(0.06) nm. 

Radiation damping implies that the initial state wave-function is not a groundstate, 
and has an average equilibrium principle quantum number n. Then our result is modified: 

Ax = VE^f (4.58) 

ESq 

n is estimated to be 

/ m 3 \ 1/2 

s ~ Ufe < 4 - 59 > 
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or, for the above parameters, n ~ 0.11 x 10 4 . Eg. ( [4.59 ) is essentially classical, and yields 
Ax ~ 2 nm, roughly consistent with the classical vertical final focus beam spot size of the 
NLC, ~ 5 nm. A more precise analysis of this latter effect is, however, certainly required. 
This may be a useful way to approach other phenomena, such as the Oide effect, in which 
large fluctuation in n in strong final focus magnets can occur, broadening the beam spot 
size. 

We have examined the quantum solution in a simple FODO synchrotron model. In 
a synchrotron information about the initial 5q is lost, and the minimal transverse beam 
spot size is: 

Ax = (4.60) 

which is 0(1) nm for most high energy synchrotrons, e.g., LEP and Tevatron, in operation 
at present. Again, a factor of n would yiedl the classical result. Presumably proton 
synchrotrons are far from the equilibrium, with n » n. 

A tantalizing question is: can quantum diffractive effects be observed? More generally, 
our discussion has been motivated by the belief that quantum optics may be the preferred 
way to analyze futuristic machines. A more general formalism more symmetrical in p± and 
x± for the study of the quantum phase space, perhaps based upon Wigner's formulation 
of quantum mechanics, is desired [15[ . 
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